Kinetics of Modulated and Ordered Structures in CuAu. 
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A continuum model derived from an atomistic Hamiltonian 
is used to examine the ordering kinetics in CuAu. A detailed 
description of the formation of the low and high tempera- 
ture ordered and modulated superlattice states is given. The 
metastability of the modulated phase at low temperatures is 
shown to severely hinder creation of the ordered superlattice. 
Formation of the modulated superlattice is shown to results 
in interesting lamellar and labyrinthine structures. 
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The vast majority of synthetic and naturally occurring 
materials contain long-lived non-equilibrium morpholo- 
gies which determine material properties. Understand- 
ing the non-equilibrium or kinetic processes that lead 
to such microscopic structures is, therefore, one of the 
most important tasks facing materials theory. The binary 
alloy CuAu is an excellent material to study such pro- 
cesses for the following reasons: the existence of a modu- 
lated superlattice at intermediate temperatures gives rise 
to a rich set of phase transformations and microstruc- 
tures; the kinetics are accessible in x-ray scattering ex- 
periments; and a continuum model for the kinetics has 
been derived from a quantum mechanical description of 
this alloy |Q. In this paper the continuum model will be 
used to provide a detailed theoretical description of the 
non-equilibrium processes that occur in CuAu. Thus pre- 
dictions that can be verified experimentally are obtained 
from a microscopic description of the alloy. 

The continuum model predicts two first order phase 
transitions; one from a disordered superlattice to a mod- 
ulated superlattice at T — Tmd and the other from the 
modulated phase to an ordered superlattice at T = Tom 
(where Tom < Tmd)- The purpose of this paper is to 
determine the transient morphologies that arise following 
instantaneous temperature quenches that bring the sys- 
tem from one equilibrium state to another. The results 
of this study indicate that: the ordered superlattice is 
difficult to form just below Tom due to the creation of a 
long-lived metastable modulated phase; the growth of an 
ordered phase far below Tom is hindered by the presence 
of 'terminal' droplets; and the nucleation of the modu- 
lated phase from the low and high temperature phases 
leads to labyrinthine and lamellar patterns respectively. 

The appearance of the modulated superlattice was pre- 
dicted from an atomistic model by Chakraborty and Xi 



Q. In this work, effective medium theory |2|^] was em- 
ployed to obtain an approximate classical Hamiltonian 
from a full quantum mechanical model. A mean field 
analysis of this Hamiltonian was used to obtain the free 
energy as a function of the average sublattice concentra- 
tion (77) along one degenerate ordering direction. This 
degeneracy can lead to interesting effects but will not be 
considered here. The free energy can be written, 



F{rj) =T a + (2Tr/a) d J dr[a T r, 2 - wf + vrf 



e|V ± r;| 2 + e|V||| 2 + /(Vir ? ) 
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where, a T = 0M2ATmeV/K, AT = T 
T = 1219K, u = 5.25meV, v = 6.1meV, e = 
25.5meV/(2Tr/a) 2 and / = 195.5meV/(2n/a) 4 and a is 
the shortest lattice parameter in the tetragonal struc- 
ture of the low temperature CuAu(I) ordered phase. The 
symbols Vj_ and V|| refer to gradients in the plane per- 
pendicular and parallel to the ordering direction respec- 
tively. At low, intermediate and high temperatures T 
is minimized by ordered (i.e., r\{x) = r\ ^ 0), mod- 
ulated (e.g., i~)(x,y,z) = r)(x + X,y,z)) and disordered 
(i.e., r)(x) = 0) states respectively. The transitions from 
the ordered to modulated phase and from the modulated 
to disordered phase are both first order and respectively 
occur at AT OM ~ —8.0K, and AT M d « 44.4K" ||,|. 
The spinodal temperatures Tq and T M are defined as the 
temperatures above which the ordered and modulated 
phases are unstable, while T5 is defined as the tempera- 
ture below which the disordered state is unstable to the 
modulated phase. The values of these temperatures are 
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T +35.7K, T M = T +52.6i4TandT| = T o +20.0K. 
The kinetics are assumed to be relaxational and driven 
by minimization of the free energy: 



d-q/dt = -TST/Srj + c, 



(2) 



where £ is a random Gaussian noise term with corre- 
lations, (((f,t)C(f,f)) = 2k B TTS(r - F)5(t - t') and 
kg is the Boltzmann constant. It is convenient to intro- 
duce the scaled variables; x = r/X and r = t/t Q where 
A = y/2f/e and t a = {a/2w) d 2f /(e 2 T) to obtain; 

dri/dr = -(7 AT + 2Vi + V| - 2V 2 ,)t7 



u'rf 
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Here, 7 = Aarf/e 2 , u' = 
{v{x,t)v{x' ,r')> = 2eS(x 



(3) 
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2k B T[a/2ix\ d [2f /e 2 }[e/2f] d / 2 . In these rescaled units the 
modulated wavelength is simply 2-7T. To simplify calcula- 
tions ordering in the parallel direction will be neglected. 

The kinetics that follow a rapid quench are influ- 
enced by the modulated phase, even when the pre- and 
post-quench temperatures are not within the modulated 
regime. For example, consider the instantaneous quench 
of a disordered state (T^^ ; > Tmd) into the ordered 
regime (T/j no ; < Tom)- The initial response will be to 
form a modulated type structure, since linear stability 
analysis predicts that the mode with q = 1 has the largest 
growth rate ||. For Tf ina i < Tom the growth of the 
q = 1 mode is much more rapid than q = 0, thus a 'mod- 
ulated' structure is quickly formed. Further growth of the 
q = mode (i.e., the ordered phase) is inhibited by the 
metastability of the modulated phase. Thus for quenches 
just below Tom a long-lived transient modulated phase 
should develop. 
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FIG. 1. In this figure the terminal droplet size is shown as 
a function of temperature. The inset depicts the velocity of 
the droplet radius as a function of R as described by Eq. (4). 

For Tfi na i "C Tom the difference between the q = 1 
and q = linear growth rates is insufficient to cre- 
ate a modulated structure. The subsequent dynamics 
do not however reduce to that observed in standard or- 
der/disorder transitions in which the average ordered do- 
main size grows as t 1 / 2 |||7| and droplets shrink at a rate 
proportional to the curvature. Here a droplet refers to 
a spherical (or circular in 2-d) regime of one phase (e.g., 
i] = —r] ) embedded in the other (e.g., t] = +r) )- In 
CuAu, large droplets shrink but do not disappear. This 
effect can be traced to the fourth order spatial derivative 
in the dynamics or equivalently, to the existence of the 
length scale describing the modulated phase, and to the 



two dimensional nature of the ordering process. 

To illustrate this effect, consider the dynamics of a 
single droplet of radius R in the limit R 3> W, where W 
is the domain wall thickness and is of the order 1 . In this 
limit 77 can be written T)(r— R(t)) ~ r] ld (r — R(t)) , where 
r) ld {x) is the one-dimensional solution of 5J-/5i] = with 
boundary conditions nix = ±00) = ±n . Substituting 
this result into Eq. (|3|), multiplying by {drfd/dr) and 
integrating from r — to r — 00 gives, 

dR/dT = 2(l + a 3 /a 1 )R+l/R 3 , (4) 

where cr, = J °° dr(dn 1 d/dr)(d' l r] 1 d/dr' 1 ). The quantity 
2(1 + 03/01) is less than zero for T < Tom- Thus the 
droplet will shrink for R > Rt and grow for R < Rt, 
where Rt is the 'terminal' droplet size and is equal to 
R T = 1/V-2(1 + 0-3/V1). The growth rate of R is il- 
lustrated in the inset of Fig. (1). While this calculation 
correctly identifies the existence of a terminal droplet, the 
accuracy is hindered by the approximation R 3> W. To 
overcome this obstacle Rt was determined numerically 
and the results are shown in Fig. (1). 
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FIG. 2. In this figure the linear dispersion for the position 
of an antiphase domain boundary is shown as a function of 
wavevector. The solid and open points represent w q in the 
absence and presence of droplets respectively and the dashed 
line corresponds to in, ~ q 2 . In the inset figures (a) and (b) 
show the interface relaxation in the absence of droplets and 
in (c) and (d) in the presence of droplets at the same times. 

The terminal droplets influence the late stage dynam- 
ics by interfering with the motion of antiphase domain 
walls or interfaces. For comparison it is useful to first 
consider the interface motion in the absence of droplets. 
Consider a slowly varying interface such that tj(x, y, t) w 
rj for x > h(y,r) and r)(x,y,T) « — r/ for x < h(y,t), 



2 



thus defining the position of the antiphase domain wall 
as h(y,t). If h(y,t) varies slowly in space it is simple to 
show {§ that dh/dr = [-2(1 + cr 3 M)V 2 + V 4 ]h using 
the methods discussed in the previous paragraph. The 
solution in Fourier space is then; h(q,r) — e~ WqT h(q, 0), 
where w q = — 2(1 + 0-3/(71 )q 2 + q 4 . In the long wavelength 
limit (q <C 1) this reduces to the standard result for or- 
der/disorder transitions (i.e., w q ~ q 2 ). To determine 
the dispersion relationship (i.e., w q ) in the presense of 
terminal droplets numerical simulations were conducted 
for the configurations shown in the inset of Fig. (2). The 
results of these calculations (see Fig. (2)) indicate that 
the relaxation of the antiphase domain walls is restrained 
by the terminal droplets. This simple effect will strongly 
alter the late stage dynamics as the droplets tend to ac- 
cumulate at interfaces Q. This accumulation is not due 
to an attraction between drops and interfaces, but rather 
by droplets getting swept up by a relaxing interface. 

The transient morphologies that emerge in the mod- 
ulated regime depend on both the pre- and post-quench 
states. If the pre-quench temperature is above Tmd and 
the post-quench temperature is between Tom and Tj, 77 
is linearly unstable and a convoluted modulated struc- 
ture will quickly emerge. Similar behavior will be ob- 
served if the pre-quench temperature is below Tom and 
the post-quench temperature is between Tq and Tmd- 
In contrast, the kinetics are dominated by nucleation and 
growth when the pre-quench state is metastable at the 
post-quench temperature. In addition growth of the nu- 
cleated droplets is not a simple process since there is 
an internal structure within the drops. Each drop will 
contain a lamellar structure that defines the modulated 
phase. In general the velocity of the droplet fronts will 
depend on the orientation of lamella with respect to the 
droplet front and on the background matrix the droplet 
is growing in. To understand this effect it is useful to es- 
timate the growth velocities in directions perpendicular 
and parallel to the lamella as the drop grows into either 
disordered or ordered backgrounds. 

First consider the invasion of the modulated phase into 
the ordered or disordered phase in a direction parallel to 
the lamella. The velocity of this front will be denoted vP 

and vP for propagation into the ordered and disordered 
phases respectively. To estimate vP , 77 can be approxi- 
mated as {ri + ri ld (y))F° (x,t) — r) , where r\ is the value 

of 77 in the ordered phase, i] ld (y) is the one dimensional 
modulated solution (i.e., rj ld (y) « A sin(g y)) and the 
overlap function FP{x, t) takes the values 1 and in the 
modulated and ordered phases respectively. Substituting 
this approximation into Eq. (^), multiplying by r] ld (y), 
and averaging over one wavelength in y gives, 



5v> v 2 (4 v 2 + <3u>A 2 )(FP) 4 
57/(77* + 3tj A 2 /2 + ^ 4 /8)(if ) 5 . 



dFP/dr = 



( 1 - 7 AT + 3rfcu' - 577*7/ 
+ 277 o 2 (107;'77 o 2 -37 J ')(F | f) 2 
+ 3(77 2 + ^ 2 /4)(7 J '- 107 7o V)(F | f) 3 



)*lf 



(5) 



if r] ld is approximated by a single Fourier mode. The 
front velocity can be obtained by expanding around the 
zero velocity limit which occurs at T = Tom- Substitut- 
ing the approximation, F® (x, t) s=s fff(x — vPt) (were ff? 

- Tom, with boundary 
(00) = 0)) into Eq. (g) 

gives, 



is the solution of Eq. (p) at T 



conditions fP{— 00 



1 and fP 



and integrating over dfP (x) / dx 



0.045(T- 



Tom)/0\\ i 



(6) 



where aP = 



J dx(dfP(x)/dx) 2 . The same tech- 
niques can be used to estimate vP and leads to vP = 
0.025(T M r> - T)/aP, where aP = J dx(dfP (x)/dx) 2 , 
and fP(x) is the solution of Eq. with r\ — at 
T = Tmd and boundary conditions fP(— 00) = 1 and 
fP (00) — 0. Numerically determining aP and aP gives, 
v° « 0.13(T - Tom) and 7;^ « 0.12(T MD - T). 

Next consider growth in the direction perpendicular to 
the lamella. For these calculations the velocities will be 
denoted vj_ and for growth into the ordered and dis- 
ordered regimes respectively. To estimate vP_, 77 can be 
approximated as; rj(x, r) « r] ld (x)Ff (x, r). Substituting 
this expression into Eq. (||), multiplying by rj ld (x), aver- 
aging over one wavelength and assuming that F 1 ?(x 1 t) 
varies slowly in space relative to the modulated wave- 
length gives, 



dFf/dr= ( 1-7AT + 49,, 
+ (37i' J 4 2 /4)(Ff ) 3 



{5v'A 4 /8){Fff 



(7) 



if a one mode approximation for 77 ld is used. Employing 
the method outlined above gives, vP^ — 0.025(Tmd ~ 
T)/aP_. Determining aP_ numerically gives m 
0.19(Tmd — T). It is much more difficult to estimate vP_ 
since the solution can not be written as an overlap func- 
tion times the modulated solution. Numerical attempts 
to determine vP_ indicate this velocity is extremely small, 
in fact no growth in this direction was observed. Thus 
to the accuracy of the current work « 0. The rea- 
son for this slow growth is that the modulated solution 
is very similar to the amplitude of the ordered solution 
and consequently the spatial transition from the ordered 
to modulated phase can be extremely sharp. In the limit 
of an infinitely sharp interface a — > 00 and v — » 0. 

The results of these calculations indicate that for 
quenches from the disordered phase the droplets will 
grow asymmetrically in a direction perpendicular to the 
lamella, leading to elongated lamella structures since 
vP < v®. In contrast, for quenches from the ordered 

phase vP > v± . Thus a small spherical nucleated droplet 
will tend to sprout arms which will invade the neighbor- 
ing territory. This growth leads to a labyrinth type struc- 
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ture similar to that observed in reacting chemical fronts 

To illustrate the remarks of the preceding paragraphs 
several numerical simulations were conducted, which are 
illustrated in Fig. (3). The numerical algorithm is dis- 
cussed in a previous paper Q and is based on the cell 
dynamics method ||. Figures (3a) and (3b) respec- 
tively show transient patterns that emerge from quenches 
from the disordered regime to just below and far below 
Tom- The highly interconnected morphology shown in 
Fig. (3a) is a long-lived metastable modulated struc- 
ture that typically evolves only near defects. The deeper 
quench (Fig. (3b)) shows that regions of high interface 
curvature coincide with large concentrations of terminal 
droplets, thus indicating the droplets inhibit interface re- 
laxation and domain growth. Figures (3c) and (3d) are 
transient patterns obtained by the nucleation of the mod- 
ulated phase from quenches from the disordered and or- 
dered regimes respectively. These figures show the lamel- 
lar and labyrinthine type structures predicted in the pre- 
ceding paragraphs. 




FIG. 3. Figures (a) and (b) are the transient patterns that 
emerge following quenches from a disordered state. Figures 
(a) and (b) are at r = 182 and r = 329 at quench tempera- 
tures of T = Tom - ±K and T = Tom - AAK respectively. 
Figure (c) is a pattern obtained by quenching from the disor- 
dered state into the modulated regime at T = Tom + S9.5K 
at t — 57.3. Figure (d) is a pattern obtained by quench- 
ing from the ordered state into the modulated regime at 
T = Tom + 38K at r = 84.0. 

In summary several predictions for the kinetics of 
phase transformations in CuAu have been made. For 
quenches from the disordered to just below the or- 
dered/modulated transition it is predicted that linear 



dispersion will produce a modulated morphology. This 
modulated structure will be long-lived since the mod- 
ulated phase is metastable in the ordered regime. For 
quenches far below Tom it was found that the formation 
of the ordered phases are strongly inhibited by the ap- 
pearance of terminal drops. Growth in this region will 
be much slower than standard order/disorder transitions 
since the droplets collect at antiphase boundaries and in- 
hibit motion. For quenches into the modulated phases it 
was found the initial state played an important role in 
determining the transient morphologies. Quenches from 
the ordered and disordered region respectively produced 
labyrinthine and lamellar type structures. 

Finally it is interesting to speculate on the asymp- 
totic dynamics in the modulated regime. In this case 
the dynamics are controlled by the elimination of defects 
and the relaxation of the orientation of the individual 
lamella. Similar dynamics occur in the Swift Hohenberg 
|Tof model of Rayleigh-Benard convection and lead to the 
growth of domains of the same orientation growing at at 
rate of i 1 / 4 in the presence of thermal fluctuations and 
t 1 / 5 in the absence of fluctuations jnj. It is likely the 
same growth exponent would be measured in this sys- 
tem. 
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